超音波影像處理

专利摘要:
提供用於處理超音波影像，特定言之不均勻組織之影像的方法、系統及電腦程式產品。方法包括：在已接收超音波影像中對組織分段；產生以組織區域中像素為中心的像素之遮罩，其中遮罩中的像素處於中心像素之給定分貝範圍之內且處於分段邊界之內；以及藉由估計遮罩中的像素之衰減將衰減分配至中心像素。已接收超音波影像可為原始影像，方法可包括產生組織區域中之每個像素之衰減的衰減圖的影像。或者，已接收超音波影像可為時間增益補償影像，方法可包括決定相對高的或低的衰減之區域。
公开号:TW201302161A
申请号:TW101119171
申请日:2012-05-29
公开日:2013-01-16
发明作者:Sivan Gieichman；Eugene Walach
申请人:Ibm；
IPC主号:G01S7-00

专利说明:
超音波影像處理
本發明係關於超音波影像處理之領域。特定言之，本發明係關於不均勻的組織之超音波影像處理。
先前技術及描述使用乳房組織診斷之情況。然而，所描述方法和系統亦可應用於其他形式的組織，特定言之不均勻類型的組織。
通常，使用X射線乳房攝影術偵測及診斷乳癌。此檢驗存在三個主要的問題：該檢驗將患者曝露於游離輻射；該檢驗對於具有緻密乳房之女性(一般而言年輕女性)不太有效；該檢驗不能在諸如活體組織切片或外科之手術期間使用。另一方面，B型掃描乳房超音波對於緻密乳房無害、有效，且易於在手術期間使用。然而，已發現作為乳癌之重要的指標之微鈣化在超音波影像中很少可見。此外，超音波影像之解釋，尤其在識別微鈣化的情況下的超音波影像之解釋在很大程度上取決於檢查醫生。
藉由對來自乳房超音波影像的聲音衰減繪圖，可得到能有助於尋找乳房腫塊及微鈣化的更加定量的標誌，從而減少對檢查醫生的依賴性。相同的方法可用於乳房瘤之特徵化(詳言之，癌性的情況與良性的情況之間的區分)。
超音波中之衰減為作為穿過成像媒介的距離之函數的超音波束之振幅的減小。考慮超音波中的衰減效應很重要，因為減少的訊號振幅可影響所產生影像的品質。藉由認知超音波束穿過媒介之衰減，可調整輸入訊號振幅以補償在所要成像深度上的任何能量損失。
根據本發明之第一態樣，提供用於處理藉由電腦化裝置執行的超音波影像之電腦實施方法，該方法包含以下步驟：在已接收超音波影像中對組織分段；產生以組織區域中之像素為中心的像素之遮罩，其中遮罩中的像素處於中心像素之給定分貝範圍之內且處於分段邊界之內；以及藉由估計遮罩中的像素之衰減將衰減分配至中心像素。
根據本發明之第二態樣，提供用於處理超音波影像的電腦程式產品，該電腦程式產品包含：包含有電腦可讀取程式代碼的電腦可讀取非暫時性儲存媒體，該電腦可讀取程式代碼包含：經配置以進行以下操作之電腦可讀取程式代碼：在已接收超音波影像中對組織分段；產生以組織區域中之像素為中心的像素之遮罩，其中遮罩中的像素處於中心像素之給定分貝範圍之內且處於分段邊界之內；以及藉由估計遮罩中的像素之衰減將衰減分配至中心像素。
根據本發明之第三態樣，提供用於超音波影像處理之系統，該系統包含：處理器；超音波資料接收器；用於在已接收超音波影像中對組織分段的分段組件；用於產生以組織區域中之像素為中心的像素之遮罩的遮罩組件，其中遮罩中的像素在中心像素的給定分貝範圍之內及在分段邊界之內；及用於藉由估計遮罩中的像素之衰減將衰減分配至中心像素的衰減估計組件。
在以下詳細說明中，闡明許多特定細節以便提供對於本發明的徹底的瞭解。然而，熟習此項技術者將瞭解，可在沒有該等特定細節的條件下實施本發明。在其他情況下，未詳細描述熟知方法、程序及組件以便不使本發明模糊。
本文中所使用的術語僅為了描述特別實施例之目的，且不欲視為本發明之限定。除非上下文另外清楚地標明，否則如本文所用的單數形式「一」、「一個」及「該」意欲同樣包括複數形式。應進一步瞭解，當術語「包含」及/或「含有」用於本說明書中時，規定所指明特徵結構、整數、步驟、操作、元件及/或組件之存在，但是並未排除一或更多其他的特徵結構、整數、步驟、操作、元件、組件及/或上述各物之群體之存在或附加。
在以下請求項中相應的結構、材料、動作及所有方法或步驟的均等物加功能元件意欲包括用於執行與明確地申請的其它主張元素相結合的功能之任何結構、材料或動作。提供本發明之描述是為了達到說明及描述之目的，而不欲為窮盡的或限於所揭示形式的本發明。在不脫離本發明之範疇及精神的情況下，許多修改及變化將對於一般技藝者是顯而易見的。選擇及描述實施例以便最佳地解釋本發明之原理及實際應用，並使得其它一般技藝者瞭解用於具有各種修改的各種實施例的適合於所考慮的特殊用途的發明。
所描述超音波影像處理在成像應用中允許用於(例如)乳房組織之不均勻組織的衰減方法。所描述超音波影像處理首先將衰減成像方法應用於不均勻組織，其次使得經處理(非原始)之超音波影像能够實施，允許該方法用於經處理影像與原始影像兩者。
為了估計組織衰減，不得不識別每個像素周圍的統一均勻組織區塊。然後計算全部區塊的平均衰減並將平均衰減分配至中心像素。以目前技術方法，此等統一區塊藉由尋找具有與中心像素之值相對類似的值的長游程像素而識別。此方法良好地用於可容易發現長均勻游程的均勻器言(諸如肝)。然而，在不均勻組織中，典型的範圍內游程頗短。因此可靠的衰減估計變得不可能。
為了緩和該問題，描述了用於包含範圍內與範圍外像素兩者的區塊之估計方法。該方法藉由局部的及全局的分段方法之組合進一步加強。
在原始超音波資料上實施所描述演算法將提供對真實組織衰減值之繪圖。該方法提供定量資料(與由標準超聲波硬體提供的定性影像相反)。不幸的是，在目前技術系統中，原始資料通常為不可用的。確實，系統供應者力圖產生均勻顯示。相應地，該等系統供應者趨於放大來自身體內部深處的組織之微弱訊號。此情況藉由應用放大(來自較深處的)遲訊號的所謂TGC(時間增益補償)曲線達成。通常，影像處理在連接至醫院PACS(圖像歸檔和通訊系統)的顯示工作站完成，在該PACS中TGCe計不可實行。然而，如下文所述，可計算每一個像素相對於像素周圍的相對衰減。因此，保留了衰減影像之診斷值。
參看第1圖，流程圖100圖示所描述超音波影像處理之實施例。
在101處接收超音波影像，及完整的超音波影像中所關注的組織之分段可執行於102處。例如，超音波影像可為忽視來自上面的脂肪組織及來自下面的胸腔的乳房組織及分段。或者，超音波影像可為肝、前列腺或其他組織。
在103處，定義組織區域中每個像素周圍的區塊。在104處，產生每個區塊之遮罩，該每個區塊指示在中心像素給定分貝範圍內之像素，此步驟的一個實施例描述於以下第4圖中。在105處，然後可根據在步驟102中識別的宏段切割遮罩。
在106處可單獨根據遮罩上之像素估計每個區塊之衰減，且在107處可將該值分配至中心像素。
在108處決定資料是否為原始資料。若資料為原始資料，在109處可產生衰減圖之影像且可識別可疑區域。
在110處，若資料不是原始資料且影像具有時間增益補償(TGC)，可只標記具有高的或低的衰減的像素。即，可丟棄具有接近零的衰減的全部像素。以此方式可找到組織中的可疑區域。
參閱第2圖，方塊圖圖示超音波影像處理系統200。影像處理系統200包括超音波資料接收器201，該超音波資料接收器201可接收原始影像資料或具有時間增益補償之影像。影像處理系統200可包括用於在完整的超音波影像中分段所關注的組織的分段組件202。
影像處理系統200亦可包括遮罩組件210。遮罩組件210可包括用於定義組織區域中每個像素周圍的區塊之區塊定義組件211。遮罩組件210可包括用於產生在中心像素之給定分貝範圍內的像素之每個區塊之遮罩的遮罩產生組件212。遮罩組件210可進一步包括遮罩切割組件213，該遮罩切割組件213用於根據在分段組件202識別的宏段切割遮罩。
影像處理系統200可包括衰減估計組件220，該衰減估計組件220用於根據遮罩上的像素估計每個區塊之衰減並將此值分配至中心像素。
影像處理系統200可包括用於產生衰減圖之影像以輸出及分析的衰減繪圖組件221。衰減繪圖組件221可在原始資料接收於影像處理系統200處時使用。
影像處理系統200可包括用於標記具有衰減的已接收影像的影像標記組件222。影像標記組件222可在已接收影像為時間增益補償影像時使用。
影像處理系統可包括用於輸出經衰減之繪圖影像或具有衰減標記的影像以分析的輸出組件230。
參閱第三圖，用於實施本發明之態樣的示例性系統包括資料處理系統300，該資料處理系統300適合於儲存及/或執行程式代碼，該程式代碼包括經由匯流排系統303直接地或間接地耦接至記憶體元件的至少一個處理器301。記憶體元件可包括在實際執行程式代碼期間採用的本端記憶體、大容量儲存器及快取記憶體，該等快取記憶體提供至少一些程式代碼的暫時儲存以便減少執行期間需從大容量儲存器中擷取代碼的次數。
記憶體元件可包括系統記憶體302，該系統記憶體302為唯讀記憶體(ROM)304及隨機存取記憶體(RAM)305之形式。基本輸入/輸出系統(BIOS)306可儲存於ROM304中。系統軟體307可儲存於包括作業系統軟體308的RAM 305中。軟體應用e310亦可儲存於RAM 305中。
系統300亦可包括諸如磁性硬碟驅動之主儲存器構件311及諸如磁碟驅動及光碟驅動之輔助儲存器構件312。驅動器及驅動器相關聯的電腦可讀取媒體提供電腦可執行指令、資料結構、程式模組及用於系統300的其他資料的非揮發性儲存。軟體應用程式可儲存於主儲存構件311及輔助儲存構件312上以及系統記憶體302上。
計算系統300可經由網路配接器316在使用邏輯連接至一或更多個遠程電腦的網路環境中操作。
輸入/輸出裝置313可直接地耦接至系統或者經由中間I/O控制器耦接至系統。使用者可經由諸如鍵盤、指向裝置或其他的輸入裝置(例如微音器、操縱桿、遊戲板、圓盤式衛星電視天線、掃描器等)的輸入裝置將命令及資訊輸入系統300中。輸出裝置可包括揚聲器、列印機等。顯示裝置314亦經由諸如視訊配接器315之介面連接至系統匯流排303。
現提供示於第1圖中之流程圖之方法的更多細節。
在一個實施例中，使用經設計以識別影像之相對統一的部分的主動等高線法完成影像中分段組織之步驟(第1圖中之步驟102)。例如，在典型的乳房影像中，所關注的區域將該區域顯示為明亮的水平條紋區域。
定義影像的每一個像素周圍的區塊之步驟(第1圖中之步驟103)可選擇在中心處具有已選擇像素的環繞矩形區塊。區塊可具有9至25個像素範圍內之寬度及25至125個像素範圍內之高度，(例如)在一個實施例中，區塊可具有17x65像素之尺寸。
參閱第4圖，流程圖400顯示第1圖之遮罩產生步驟104之實施例。在401處，接收中心像素周圍的超音波影像之區塊。在402處，產生初始遮罩，該初始遮罩指示處於中心像素的給定分貝範圍之內的像素，(例如)在距離中心像素2至6分貝的範圍之內的像素。在一個實施例中，可產生在中心像素3分貝內像素之遮罩。在403處，可在遮罩中填入小的垂直空隙。
接著，在404處，可將遮罩劃分為連接組件。在405處，可移除除了中心連接組件之外的全部連接組件。在406處，可從遮罩中移除在步驟403中加入的用以填充空隙的像素。在407處，亦可從遮罩中移除具有少數非零元件的行。然後在408處輸出遮罩。
或者，不僅僅保留中心連接組件，而是分別估計每個連接組件的衰減，且最後的衰減值為全部的連接組件之平均值。
以下代碼為用於遮罩產生步驟之Matlab(矩陣實驗室)(MATLAB為MathWorks公司之商標)數值計算環境之實例：function mask=CreatMask(E,x,y)h=size(E,1)；pc=E(y,x)；中心像素mask=(E<(1.3*pc))&(E>(0.7*pc))；像素範圍mask1=mask；填充孔[labels,num]=bwlabel(~mask)；劃分為連接組件for L=1：num；查看全部連接組件Lmask=(labels==L)；當前連接組件row=find(sum(Lmask,2))；連接組件之列索引if(min(row)>1)&&(max(row)<h)&&(sum(sum(Lmask)>4)<3)mask=mask+Lmask；end end mask=bwselect(mask,x,y)；僅保留中心CC於遮罩上mask=mask.*mask1；在第1圖中之步驟105中，越過識別於第1圖中之分段步驟102的宏段之邊界的部分遮罩被移除。
對於估計每個區塊之衰減之步驟(第1圖中之步驟106)，假定用於反向散射回波振幅的簡單指數模型。在假定遮罩上的全部像素中的衰減是統一的情況下，使用擴展的普諾尼(Prony)方法之變化完成對於衰減之估計。普諾尼的方法從統一取樣的訊號中提取有價值的資訊並構建一系列阻尼複指數或正弦曲線。
P及K分別表示區塊寬度及高度。亦以Ω _i 表示遮罩之第i行，及以k _i 表示Ω _i 中非零元件之數目。在切割遮罩之步驟中，以下價值函數被最小化： 其中E _i,j 為回波振幅(超音波影像中第i,j像素之強度)；E _i 為區塊之第i行頂端之振幅；以及a為區塊之衰減。
將相對於E _i 之導數設為零得到：
使用指數之泰勒級數展開得到：=a _I -2b _i a，其中： 代入價值函數c並將相對於a之導數設為零意謂： 其中：D _i =(a _i -E _i,1,…,a _i -E _i,k ) ^T ,B _i =2．diag(Ωi)．(b _i ,b _i +a _i ,…,b _i +(k _i -1)a _i ) ^T .
在替代性實施例中，遮罩中之每行分成連續的非零段(同時忽略小孔)，且上述的衰減(α)之計算藉由對來自全部的行(增加P)的全部的分開的段求和執行。在此情況下，不假定衰減沿著遮罩是統一的，而是假定衰減為分段的常數，亦即，衰減在連續的非零段是統一的。
以下為衰減計算之Matlab實施之實例：載入影像im_name='U5cropped.tif'；im=imread(im_name)；if size(size(im),2)==3 im=rgb2gray(im)；end dim_im=size(im)；figure(1)imshow(im,[])im=double(im)；im=im/max(im(：))；方框參數w=17；h=65；w2=floor(w/2)；h2=floor(h/2)；衰減計算alpha_map=zeros(size(im))；wb=waitbar(0,['Processing image：',im_name,'...'])；產生進度條for x=1：dim_im(2)left=max(1,x-w2)；right=min(size(im,2),x+w2)；for y=1：dim_im(1)top=max(1,y-h2)；bottom=min(size(im,1),y+h2)；E=double(im(top：bottom,left：right))；dim=size(E)；mask=CreatMask(E,x-left+1,y-top+1)；if max(sum(mask))<=10 continue；end E=E.*mask；K=sum(mask)；E(：,k<4)=[]；mask(：,k<4)=[]；k(k<4)=[]；dim=size(E)；a=sum(E)./k；A=repmat(a,dim(1),1)；D=A-E；j=(0：(dim(1)-1))’；J=repmat(j,1,dim(2))；s=2*sum(J.*mask)./k；b=sum(J.*E)./k-s.*a；B_mat=repmat(b,dim(1),1)；B=2*(B_mat+J.*A).*mask；alpha_map(y,x)=trace(D’*B)/trace(B’*B)；end waitbar(x/dim_im(2))更新進度條end close(wb)關閉進度條
可使用其他替代方法估計衰減。舉例而言，可使用最大似然量法代替普諾尼方法。然而，以上實施具有優點，因為該實施不依靠關於統計資料分佈的任何基本假設。
在處理原始資料時，估計衰減之後，可產生衰減圖之影像且可找到可疑的區域(例如乳房囊腫、瘤腫塊或鈣化)。然而，在處理具有TGC之影像時，衰減估計之後的衰減圖幾乎全部為零。因此，僅僅在具有相對高的或低的表減的區域標記初始影像而不是繪製該圖。
參閱第5A及5B圖，第5A及5B圖顯示示例性影像。在第5A圖中，所關注的組織之分段圖示於500處。在此實例中，乳房組織501從以上脂肪組織502及以下胸腔區域503分段。
在第5B圖中，顯示圖示了已標記衰減之示例性影像550。在該影像550中，顯示標記於影像550上之具有高衰減552之區域的瘤551。
如熟習此項技藝者所將瞭解，可將本發明之態樣實施為系統、方法或電腦程式產品。因此，本發明之態樣可採取完全硬體實施例之形式、完全軟體實施例之形式(包括韌體、常駐軟體、微代碼等)或在本文中可全部大體上稱為「電路」、「模組」或「系統」的組合軟體及硬體態樣之實施例。此外，本發明之態樣可採取實施於具有電腦可讀取程式代碼的一或更多個電腦可讀取媒體中的電腦程式產品之形式。
可利用一或更多個電腦可讀取媒體之任何組合。電腦可讀取媒介可為電腦可讀取訊號媒體或電腦可讀取儲存媒體。電腦可讀取儲存媒體可為(例如)但不限於：電子的、磁的、光的、電磁的、紅外的或半導體系統、設備、或裝置或以上之任何適當的組合。電腦可讀取儲存媒體之更特定實例(非窮舉清單)將包括以下：具有一或更多個電線的電連接、便攜電腦磁碟、硬碟、隨機存取記憶體(RAM)、唯讀記憶體(ROM)、可擦可程式唯讀記憶體(EPROM或閃速記憶體)、光纖、便攜光碟唯讀記憶體(CD-ROM)、光儲存裝置、磁儲存器裝置或以上之任何適當的組合。在此文件之上下文中，電腦可讀取儲存媒體可為包含或儲存供指令執行系統、設備或裝置使用的或結合指令執行系統、設備或裝置使用的程式之任何有形媒體。
電腦可讀取訊號媒體可包括(例如)在基帶中或作為載波之部分包含電腦可讀取程式代碼的傳播資料訊號。此類傳播信號可採取各種形式包括(但不限於)電磁的、光的或上述任何適當的組合。電腦可讀取訊號媒體可為任何非電腦可讀取儲存媒體之電腦可讀取媒體，且該電腦可讀取訊號媒體可通訊、傳播、或輸送供指令執行系統、設備或裝置使用的或結合指令執行系統、設備或裝置使用的程式。
包含於電腦可讀取媒介上之程式代碼可使用任何適當的媒體傳送，適當的媒體包括(但不限於)無線、纜線、光纖電纜、RF等或以上之任何適當的組合。
用於實施本發明之態樣的操作的電腦程式代碼可以一或更多程式化語言之任何組合編寫，該一或更多程式化語言包括面向對象程式化語言(諸如Java、Smalltalk、C++等)及習知程序程式化語言(諸如“C”程式化語言或類似的程式化語言)。程式代碼可完全地在使用者的電腦上執行、作為獨立套裝軟體部分地在使用者電腦上執行、部分地在使用者電腦上、部分地在遠端電腦上或完全地在遠端電腦或伺服器上執行。在後一情形中，遠端電腦可經由任何類型之網路連接至使用者的電腦，任何類型之網路包括局域網(LAN)或寬域網(WAN)，或可連接至外部的電腦(例如經由使用網際網路服務提供商的網際網路)。
參照根據本發明之實施例之方法、設備(系統)及電腦程式產品之流程圖說明及/或方塊圖描述本發明之態樣。將瞭解，流程圖說明及/或方塊圖之每個方塊及流程圖說明及/或方塊圖中方塊之組合可由電腦程式指令實施。該等電腦程式指令可提供至通用電腦、專用電腦或其他可程式化資料處理設備之處理器以生產機器，以使得經由電腦或其他的可程式化資料處理設備之處理器執行的指令產生用於實施流程圖及/或方塊圖或方塊中規定的功能/動作之構件。
該等電腦程式指令亦可儲存於可引導電腦、其他可程式化資料處理設備或其他裝置以特別方式工作的電腦可讀取媒體中，以使得儲存於電腦可讀取媒體中的指令產生包括實施規定於流程圖及/或方塊圖或方塊中之功能/動作的指令之製品。
該等電腦程式指令亦可載入電腦、其他可程式化資料處理設備或其他裝置上以致使一系列操作步驟執行於電腦、其他可程式化設備或其他裝置上，產生電腦實施過程，以使得執行於電腦或其他可程式化設備上之指令提供用於實施規定於流程圖及/或方塊圖或方塊中之功能/動作的過程。
圖中之流程圖及方塊圖圖示根據本發明的各種實施例之系統、方法及電腦程式產品之可能的實施之結構、功能性及操作。在此點上，流程圖或方塊圖中之每個方塊可代表模組、部分或代碼之部分，模組、部分或代碼之部分包含用於實施規定的一或多個邏輯函數之一或更多個可執行指令。亦應注意，在一些替代性實施例中，如方塊中表明之功能可能不以圖表中表明之次序發生。例如，事實上，兩個相繼示出的方塊可實質上同時執行，或者根據所包含的功能性，有時方塊可以相反次序執行。亦應注意，方塊圖及/或流程圖說明中之每個方塊及方塊圖及/或流程圖說明中的方塊之組合可由執行所規定功能或動作或專用硬體與電腦指令之組合的基於專用硬體的系統實施。
100‧‧‧流程圖
101‧‧‧方塊
102‧‧‧方塊
103‧‧‧方塊
104‧‧‧方塊
105‧‧‧方塊
106‧‧‧方塊
107‧‧‧方塊
108‧‧‧方塊
109‧‧‧方塊
110‧‧‧方塊
200‧‧‧超音波影像處理系統
201‧‧‧超音波資料接收器
202‧‧‧分段組件
210‧‧‧遮罩組件
211‧‧‧區塊定義組件
212‧‧‧遮罩產生組件
213‧‧‧遮罩切割組件
220‧‧‧影像處理系統
221‧‧‧衰減繪圖組件
222‧‧‧影像標記組件
230‧‧‧輸出組件
300‧‧‧資料處理系統
301‧‧‧處理器
302‧‧‧系統記憶體
303‧‧‧系統匯流排
304‧‧‧唯讀記憶體
305‧‧‧隨機存取記憶體
306‧‧‧基本輸入/輸出系統
307‧‧‧系統軟體
308‧‧‧作業系統軟體
310‧‧‧軟體應用
311‧‧‧主儲存器構件
312‧‧‧輔助儲存器構件
313‧‧‧輸入/輸出裝置
314‧‧‧顯示裝置
315‧‧‧視訊配接器
316‧‧‧網路配接器
400‧‧‧流程圖
401‧‧‧方塊
402‧‧‧方塊
403‧‧‧方塊
404‧‧‧方塊
405‧‧‧方塊
406‧‧‧方塊
407‧‧‧方塊
408‧‧‧方塊
500‧‧‧方塊
501‧‧‧乳房組織
502‧‧‧脂肪組織
503‧‧‧胸腔區域
550‧‧‧示例性影像
551‧‧‧瘤
552‧‧‧高衰減
在本說明書的結束部分特別指出並清楚地主張被認為是本發明之標的。在閱讀附圖時藉由參閱以上【實施方式】可最佳地理解關於本發明之組織與操作方法兩者，以及本發明之目標、特徵及優點，在附圖中：第1圖為根據本發明之方法之流程圖；第2圖為根據本發明之系統之方塊圖；第3圖為可實施本發明之電腦系統之方塊圖；第4圖為根據本發明之方法的態樣之流程圖；以及第5A及5B圖為圖示根據本發明的成像結果之圖表。
應瞭解，為了說明之簡單與清晰起見，示於圖中之元件不必按比例繪製。例如，為清晰起見，一些元件之尺寸可相對於其他元件誇示。此外，在被認為合適的情況下，元件符號可在圖中重複以指示相應的或類似的特徵結構。
100‧‧‧流程圖
101‧‧‧方塊
102‧‧‧方塊
103‧‧‧方塊
104‧‧‧方塊
105‧‧‧方塊
106‧‧‧方塊
107‧‧‧方塊
108‧‧‧方塊
109‧‧‧方塊
110‧‧‧方塊

权利要求:
Claims (15)
[1] 一種用於處理藉由一電腦化裝置執行的一超音波影像之電腦實施方法，該電腦實施方法包含以下步驟：在一已接收的超音波影像中對組織分段；產生以一組織區域中之一像素為中心之像素之一遮罩，其中該遮罩中之該等像素處於該中心像素之一給定分貝範圍之內且處於一分段邊界之內；以及藉由e該遮罩中之該等像素之一衰減將該衰減分配至該中心像素。
[2] 如請求項1所述之方法，其中一已接收超音波影像為一原始影像且該方法包括以下步驟：產生一組織區域中之每個像素之衰減的一衰減圖之一影像。
[3] 如請求項1所述之方法，其中一已接收超音波影像為一時間增益補償影像且該方法包括以下步驟：決定相對高的或低的衰減之區域。
[4] 如請求項3所述之方法，該方法包括以下步驟：用衰減資訊標記該已接收超音波影像。
[5] 如請求項1所述之方法，其中對於一組織區域中的每個像素執行產生一遮罩之步驟。
[6] 如請求項1所述之方法，該方法包括以下步驟：定義一中心像素周圍之像素之一區塊並產生該定義區塊內的像素之一遮罩。
[7] 如請求項1所述之方法，其中產生像素之一遮罩之步驟包括以下步驟：將該遮罩劃分為連接組件並選擇連接組件之一中心區域。
[8] 如請求項1所述之方法，其中估計該遮罩中之該等像素之該衰減之步驟假定遮罩上之全部像素中之該衰減為統一的；或其中估計該遮罩中之該等像素之該衰減之步驟假定該衰減為連續的非零段上之分段常數。
[9] 如請求項1所述之方法，其中估計該遮罩中之該等像素之該衰減使用一擴展的普諾尼方法。
[10] 一種用於處理一超音波影像之一電腦程式產品，該電腦程式產品包含：一電腦可讀取非暫時性儲存媒體，該電腦可讀取非暫時性儲存媒體包含有電腦可讀取程式代碼，該電腦可讀取程式代碼包含：經配置以執行以下操作的電腦可讀取程式代碼：在一已接收的超音波影像中對組織分段；產生以一組織區域中之一像素為中心之像素之一遮罩，其中該遮罩中之該等像素處於該中心像素之一給定分貝範圍之內且處於一分段邊界之內；以及藉由估計該遮罩中之該等像素之一衰減將該衰減分配至該中心像素。
[11] 一種用於超音波影像處理之系統，該系統包含：一處理器；一超音波資料接收器；一分段組件，該分段組件用於一已接收超音波影像中的組織之分段；一遮罩組件，該遮罩組件用於產生以一組織區域中之一像素為中心的像素之一遮罩，其中該遮罩中之該等像素處於該中心像素之一給定分貝範圍之內且處於一分段邊界之內；以及一衰減估計組件，該衰減估計組件用於藉由估計該遮罩中之該等像素之一衰減將該衰減分配至該中心像素。
[12] 如請求項11所述之系統，其中一已接收超音波影像為一原始影像且該系統包括：一衰減繪圖組件，該衰減繪圖組件用於產生一組織區域中之每個像素之衰減的一衰減圖之一影像；或其中一已接收超音波影像為一時間增益補償影像且該系統包括：一影像標記組件，該影像標記組件用於決定相對高的或低的衰減之區域並用衰減資訊標記該已接收的超音波影像。
[13] 如請求項11所述之系統，其中該遮罩組件包括：一區塊定義組件，該區塊定義組件用於定義一中心像素周圍之像素之一區塊並產生該定義區塊內之像素之一遮罩；或其中該遮罩組件包括：一遮罩產生組件，該遮罩產生組件用於產生像素之一遮罩，其中該遮罩中之該等像素位於一給定分貝範圍之內。
[14] 如請求項13所述之系統，其中該標記產生組件係用於將該遮罩劃分為連接組件並選擇連接組件之一中心區域。
[15] 如請求項11所述之系統，其中該遮罩組件包括：一遮罩切割組件，該遮罩切割組件用於根據該組織之該分段切割一遮罩。

类似技术:

---

公开号 | 公开日 | 专利标题

---

TWI526196B|2016-03-21|用於處理超音波影像處理之方法、電腦程式產品與系統

---

US20200085412A1|2020-03-19|System and method for using medical image fusion

---

US9375195B2|2016-06-28|System and method for real-time ultrasound guided prostate needle biopsy based on biomechanical model of the prostate from magnetic resonance imaging data

---

JP6462340B2|2019-01-30|超音波診断装置、画像処理装置及び画像処理方法

---

CN102028500B|2014-02-19|超声波诊断装置、超声波图像处理装置、以及超声波图像处理方法

---

CN106725614B|2021-03-30|超声成像设备、超声成像方法及装置

---

JP2011224362A|2011-11-10|超音波データにおいて関心領域を決定するための方法及びシステム

---

JP2013005876A|2013-01-10|超音波診断装置

---

KR20150061989A|2015-06-05|탄성 영상 제공 방법 및 이를 위한 초음파 장치

---

JPWO2009060751A1|2011-03-24|超音波診断装置とその作動方法及び超音波画像診断プログラム

---

JP2017086896A|2017-05-25|医用画像処理装置及び超音波診断装置

---

KR20160094766A|2016-08-10|의료 영상을 디스플레이 하기 위한 방법 및 장치

---

JP2017064183A|2017-04-06|医用画像診断装置、医用画像処理装置および医用情報の表示制御方法

---

RU2636262C2|2017-11-21|Способ и система для обработки данных ультразвуковой визуализации

---

Lagendijk et al.2018|Breast and tumour volume measurements in breast cancer patients using 3-D automated breast volume scanner images

---

US10206661B2|2019-02-19|Ultrasound with augmented visualization

---

Sparks et al.2015|Multiattribute probabilistic prostate elastic registration |: Application to fusion of ultrasound and magnetic resonance imaging

---

von Lavante et al.2008|Segmentation of breast cancer masses in ultrasound using radio-frequency signal derived parameters and strain estimates

---

US10881379B2|2021-01-05|Method of visualizing a sequence of ultrasound images, computer program product and ultrasound system

---

KR101107487B1|2012-01-19|영역 변화 그래프 제공 방법 및 그를 위한 초음파 시스템

---

US20190244352A1|2019-08-08|Analyzing apparatus and analyzing method

---

KR20170086311A|2017-07-26|의료 영상 장치 및 그 동작방법

---

JP2012125297A|2012-07-05|超音波診断装置、画像処理装置及びプログラム。

---

US9999402B2|2018-06-19|Automatic image segmentation

---

JP2017038705A|2017-02-23|超音波診断装置、及び超音波診断装置の制御方法

同族专利:

---

公开号 | 公开日

---

DE112012002337T5|2014-03-06|

---

TWI526196B|2016-03-21|

---

GB2504631B|2017-07-19|

---

WO2012163613A1|2012-12-06|

---

US20120308103A1|2012-12-06|

---

US8837792B2|2014-09-16|

---

CN103649771B|2016-08-17|

---

DE112012002337B4|2017-07-06|

---

CN103649771A|2014-03-19|

---

GB2504631A|2014-02-05|

---

GB201318423D0|2013-12-04|

引用文献:

---

公开号 | 申请日 | 公开日 | 申请人 | 专利标题

---

TWI485420B|2013-09-27|2015-05-21|Univ Nat Taiwan|超音波影像補償方法|IL96842D0|1990-12-31|1991-09-16|Ibm Israel|Image processing system for producing an attenuation map of a scanned image|

---

US6117080A|1997-06-04|2000-09-12|Atl Ultrasound|Ultrasonic imaging apparatus and method for breast cancer diagnosis with the use of volume rendering|

---

US8834371B2|2005-03-31|2014-09-16|Kabushiki Kaisha Toshiba|Ultrasound diagnostic apparatus and ultrasound image processing program|

---

US20080031506A1|2006-08-07|2008-02-07|Anuradha Agatheeswaran|Texture analysis for mammography computer aided diagnosis|

---

US8824762B2|2010-10-22|2014-09-02|The Johns Hopkins University|Method and system for processing ultrasound data|US9999402B2|2014-07-21|2018-06-19|International Business Machines Corporation|Automatic image segmentation|

---

JP2017063926A|2015-09-29|2017-04-06|セイコーエプソン株式会社|超音波受信信号補正装置、超音波測定装置及び超音波受信信号補正方法|

---

KR20170047873A|2015-10-26|2017-05-08|삼성메디슨 주식회사|초음파 영상 장치 및 그 제어 방법|

---

JP6705261B2|2016-04-01|2020-06-03|セイコーエプソン株式会社|画像生成装置および画像生成方法|

法律状态:

优先权:

---

申请号 | 申请日 | 专利标题

---

US13/151,303|US8837792B2|2011-06-02|2011-06-02|Ultrasound image processing|

---